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The distribution of FCoV Types I and II in a Portuguese cat population was studied by a RT- 
PCR assay targeting the 3’-end of the viral RNA. For a period of 3 years, 120 samples were 
collected and 57 were found positive for FCoV RNA. Within the positive samples the 
presence of FCoV Type I was found in 79%. Type II was only detected in 3.5% in animals 
with Feline Infectious Peritonitis. The remaining 17.5% could not be differentiated. These 
viral sequences, comprising a region within gene S were further subjected to a 
heteroduplex mobility assay (HMA) detecting the presence of viral quasispecies in 17% of 
the samples. Phylogenetic analysis for FCoV Type I revealed high genetic diversity 
between the Portuguese sequences and other previously characterized strains, while 
Type II tree showed a higher genetic homogeneity. This study confirmed the presence of 
FCoV Types I and II circulating in Portugal and detected high genetic diversity between 
circulating strains suggesting that the virus persists within the host as mixed viral 
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1. Introduction 


Feline Coronaviruses (FCoVs) are members of genus 
Coronavirus, Family Coronaviridae. Two variants of FCoV 
have been described, Feline Enteric Coronavirus responsible 
for an enteric disease and Feline Infectious Peritonitis 
Virus (FIPV) the etiological agent of FIP, a fatal disease of 
domestic and wild cats. FIP development in a small 
percentage of animals infected with FCoV is related to the 
occurrence of genomic mutations supported by the close 
genetic relationship between the higher virulent virus and 
the parental FCoV (Herrewegh et al., 1997; Vennema et al., 
1998; Vennema, 1999). Additionally, FCoVs were classified 
as Types I and II according to their serological and genomic 
features. Type I is strictly feline. Type II FCoV, although 
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sharing with Type I a high degree of similarity in the 3’-end 
of the genome, is antigenic and genetically more related to 
Canine Coronavirus (CCoV) (Hohdatsu et al. 1991), 
suggesting that Type II resulted from a homologous and 
independent recombination event within the S and the M 
genes of Type I FCoV and CCoV (Herrewegh et al., 1998). 
The current information on the geographical distribution 
of FCoV indicates a high incidence of Type I in Europe, 
Japan and USA (Addie et al., 2003; Hohdatsu et al., 1992). 
However in Japan and Taiwan detection of FIP was 
associated with Type II FCoV (Hohdatsu et al., 1992; Lin 
et al., 2009). Presently there is no available information 
about the molecular epidemiology of FCoV in Portugal. Due 
to the pathogenic impact of this virus in the feline health it 
was our purpose to investigate the presence of FCoV Types 
in the Portuguese cat population using a RT-PCR-nested 
PCR assay previously reported (Addie et al., 2003) and to 
assess its genetic diversity by heteroduplex mobility assay 
(HMA) and phylogenetic analysis. 
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2. Materials and methods 
2.1. Sampling collection 


The biological materials included rectal swabs, 
stool, plasma, thoracic and ascite fluids collected from 
animals attending the Veterinary Hospital of the Faculty 
of Veterinary Medicine, Lisbon. 


2.2. Detection of FCoV RNA 


FCoV RNA was extracted using the QIAmp Viral RNA or 
Qiagen RNeasy mini kits (Qiagen), following the manu- 
facturer’s instructions. A fragment of 177 bp located in the 
3’-UTR was amplified by RT-PCR followed by a nested PCR 
using two primer pairs described by Herrewegh et al. 
(1995). The 25 tL reactions were amplified with one step 
RT-PCR Fidelity Taq 2x Master Mix (GE Healthsciences), 
50 pmol of each primer and 10 wL of RNA sample. The 
nPCR was performed with PCR Fidelity Taq 2x Master Mix 
(GE Healhsciences), 50 pmol of each primer and 5 wLof the 
previous amplification in a 25 wL reaction 


2.3. Differentiation between Types I and II FCoV 


Positive FCoV RNA samples were subjected to a RT-PCR 
assay to amplify a 376-283 bp fragment located in the 3’- 
region of the spike gene (S), followed by a nPCR to amplify a 
360-218 bp internal fragment, for differentiation of Type I 
from Type II FCoV according to the fragment length (Addie 
et al., 2003). The amplification was performed with one 
step RT-PCR Fidelity Taq 2x Master Mix (GE Healths- 
ciences), 5 wL of RNA sample and primers Iffs, Icfs, lubs ina 
25 wL reaction, followed by a nested PCR with PCR Fidelity 
Taq 2x Master Mix (GE Healthsciences), primers nlffles, 
nicfs, nlubs (Addie et al., 2003) and 2.5 wL of the previous 
reaction in 25 wL final volume. The amplicons were 
analyzed in a 2.5% agarose gel stained with 0.5 g/mL of 
ethidium bromide in 1~x TBE. 


2.4. Heteroduplex mobility assay (HMA) and cloning 


The amplicons were subjected to a HMA before cloning 
(Delwart et al., 1993). Briefly, 10 wL of the second round S 
amplicons were combined with 1 wl of 10x annealing 
buffer (1M NaCl, 100mM Tris pH 7.8, 20mM EDTA), 
heated to 94°C/5 min, cooled on ice, added of 1 pL of 
loading buffer (50% glycerol, 0.01% of xyleno cyanol) and 
loaded in a 5% polyacrylamide gel (29.2:0.8 acrylamide:- 
bisacrylamide) poured into plates of 15cm x 15cm. The 
electrophoresis was carried at a constant 150 V/3h ina 
vertical gel apparatus (Biorad) in 1x TBE. After electro- 
phoresis the gel was stained with 0.5 wg/mL of ethidium 
bromide in 1x TBE and observed in a UV transiluminator. 
The visualized homoduplexes bands were cuted, eluted at 
56 °C/5 min in 200 «tL H20 of which 1 tL was re-amplified 
using the nPCR primers set described above. The resulting 
amplicons were cloned in pGEM®-T Easy Vector System 
(Promega) according to the manufacture’s conditions and 
sequenced by STAB Vida Lda, Portugal. The nucleotide 
sequences were submitted to Genbank and were assigned 
the accession numbers EU327694-EU327752. Sample 


identification was included in Table 1. FCoV. Types I and 
II assignment based on the amplicons length were in 
accordance with the nucleotide sequence of each sample. 


2.5. Phylogenetic analysis 


The nucleotide sequences were aligned (PILEUP from 
GCG Wisconsin Package v. 10.0), including sequences of 
FCoV Type I (KU_2 AN-D32044; UCD1 AN-AB088222; Black 
AN-AB088223; C2_0301UK AN-AY159737; J8_0196UK AN- 
AY159758; FQ1_1099UK AN-AY159754; H9_1100UK AN- 
AY159770; FUI_1191UK AN-AY159735) and FCoV Type II 
(79_1683 AN-X80799; 79_1146 AN-X06170; G6_1095 AN- 
AY159784). Due to its recombinant status FCoV Type II 
sequences were aligned and submitted to phylogenetic 
analysis separately from Type I sequences. 

The phylogenetic relationship between sequences was 
inferred by genetic distances between pairs of sequences 
using Kimura’s two parameter method (DNADIST, Phylip 
Package). The values of transition/transversion ratio and 
nucleotide frequencies were calculated from the data sets 
(TREE-PUZZLE) (Strimmer and von Haeseler, 1997). The 
resulting trees were constructed by Neighbor Joining 
(NEIGHBOR, Phylip Package) and the branching order 
reliability was evaluated by bootstrap analysis of 1000 
replicates (SEQBOOT, Phylip Package). To maximize 
the phylogenetic information of the nucleotide align- 
ment closely related sequences were excluded but were 
indicated between brackets following the sequence 
chosen to represent each set of viruses as summarized; 
65LVI_05 [74LVI_05, 81LVI_05, 86LVI_05, 87LVI_05]; 
14LVI_05 [10LVI_O5a]; 19LVI_04a [19LVI_04b]; 
58LVI_05 [59LVI_05a/b]; 1Lg_06 [1Fg_O6a]. 


3. Results 
3.1. Detection of Types I and II FCoV 


For 3 years 120 samples were collected from domestic 
cats to detect FCoV RNA due to clinical suspicion of FIP and 
to confirm viral excretion from healthy animals. Fifty- 
seven samples were found FCoV positive (47.5%) and 
afterwards processed for Types I and II differentiation. Out 
of those 57 positive samples, 45 (79%) were classified as 
Type I FCoV and were obtained from healthy animals or 
with FIP. Only two samples (3.5%) from animals with FIP, 
were classified as Type II. No amplification was obtained 
from ten samples (17.5%). 


3.2. Detection of viral quasispecies 


The presence of viral quasispecies within the animal, 
was evaluated by HMA for each sample (Fig. 1). This assay 
detects the presence of fully (homoduplexes) and partially 
(heteroduplexes) complementary DNA chains formed 
between highly related sequences during the process of 
denaturing and annealing. Homoduplexes molecules tend 
to migrate faster, while heteroduplexes show a reduced 
mobility in a non-denaturing PAGE. The observation of two 
or more homoduplexes bands, plus the corresponding 
heteroduplexes bands indicates the presence of different 


Table 1 
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Identification and origin of the Portuguese FeCoV sequences, amplified from healthy and FIP (S) animals. The genetic diversity of the amplicons was 
evaluated by heteroduplex mobility assay and each positive sample was re-amplified. The number of corresponding sequences is indicated (a-c). 


Identification Heteroduplexes Origin 

2004 

13LVI_2004 Positive (a;b) Lisbon, Portugal 
19LVI_2004 Positive (a;b) Lisbon, Portugal 
20LVI_2004 Negative Lisbon, Portugal 
22LVI_2004 Positive Palmela, Portugal Sharing household 
23LVI_2004 Negative Palmela, Portugal 

25LVI_2004 Positive (a;b) Palmela, Portugal 

2005 

3LVI_2005 Negative Lisbon, Portugal 
1OLVI_2005 Positive (a;b) Lisbon, Portugal Sharing household 
11LVI_2005 Positive (a;b) Lisbon, Portugal 

14LVI_2005 Negative Lisbon, Portugal 

27LVI_2005 Negative Lisbon, Portugal 
33LVI_2005 Positive (a;b) Rio Mouro, Portugal Multi-cat environment 
37LVI_2005 Negative Lisbon, Portugal 
38LVI_2005 Negative Lisbon, Portugal 
A8LVI_2005 Negative Lisbon, Portugal 
50LVI_2005 Negative Lisbon, Portugal 
58LVI_2005 Negative Lisbon, Portugal 
59LVI_2005 Negative Lisbon, Portugal 
71LVI_2005 Negative Lisbon, Portugal 
S8OLVI_2005 Negative Vila Nova de Gaia, Portugal 
65LVI_2005 (S) Negative Lisbon, Portugal Sharing household 
7ALVI_2005 Negative Lisbon, Portugal 

81LVI_2005 Negative Lisbon, Portugal 

86LVI_2005 Negative Lisbon, Portugal 

87LVI_2005 Negative Lisbon, Portugal 

88LVI_2005 Negative Lisbon, Portugal 

9OLVI_2005 Negative Lisbon, Portugal 

93LVI_2005 Negative Lisbon, Portugal 

9ALVI_2005 Negative Lisbon, Portugal 

96LVI_2005 Negative Lisbon, Portugal 

99LVI_2005 Positive (a;b) Lisbon, Portugal 
102LVI_2005 Negative Lisbon, Portugal 
2006 

12LVI_2006 Negative Lisbon, Portugal 
13LVI_2006 Negative Lisbon, Portugal 
15LVI_2006 Negative Lisbon, Portugal 
19LVI_2006 Negative Lisbon, Portugal 
20LVI_2006L Negative Lisbon, Portugal 
20LVI_2006PBMC Negative Lisbon, Portugal 
20LVI_2006Z Negative Lisbon, Portugal 
24LVI_2006 Negative Lisbon, Portugal 
24LVI_2006D1 Negative Lisbon, Portugal 
1FG_2006 (S) Positive (a;b) Lisbon, Portugal 
1LG_2006 Negative Lisbon, Portugal 
28LVI_2006 Negative Lisbon, Portugal 
29LVI_2006 Negative Lisbon, Portugal 
32LVI_2006 (S) Positive (a;b,c) Lisbon, Portugal 
38LVI_2006 Negative Lisbon, Portugal 


DNA populations within the sample, named viral quasis- 
pecies and mirror the in vivo genetic diversity of the virus 
(Delwart et al., 1993). 

Within the 47 S gene amplicons 37 samples (78.7%) 
showed no genetic diversity whereas in 10 samples (21.3%) 
two homoduplexes bands were observed and individually 
amplified. Comparison of the nucleotide sequences of the 
heteroduplexes positive samples revealed similarities of 
90-99% along the 320nt in Type! and 170nt in Type II FCoV. 
Although very similar, six sequences (13LVI_04b; 
19OLVI_04b; 10LVI_05a; 11LV_05a; 33LVI_O5b; 99LVI_O5b) 
showed a 5/-end variation length of 28-56 nucleotides 
justifying the homo/heteroduplexes formation. Sequences 


33LVI_05a/b, 25LVI_04a/b, 10LVI_05a/b and 11LVI_05a/b 
were amplified from animals with multi-cat environment 
(Table 2). 


3.3. Phylogenetic analysis 


The genetic relationships between the Portuguese 
sequences of FCoV and previously assigned sequences 
from isolates of Type revealed a heterogeneous cluster- 
ing of sequences with several tree branches supported 
by high bootstrap values (Fig. 2). Sequences from 
distinct viral quasispecies obtained from the same 
animal as well as sequences from cats sharing the same 
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Fig. 1. Heteroduplex mobility assay of S amplicons from samples 
1OLVI_05 (1); 11LVI_O5 (2); 3LVI_O5 (3); 25LVI_04 (4); 23LVI_04 (5); 
22LVI_04 (6); ZOLVI_04 (7); 10LVI_04 (8); 13LVI_04 (9); 12LVI_04 (10), 
included in Type I FCoV. Lanes 1, 2, 4, 6, 8 and 9 (+) show the presence of 
two homoduplexes bands (—), corresponding to fully complementary 
DNA chains. The heteroduplexes bands (*) consistent to partially 
complementary DNA chains with reduced mobility are clearly visible 
in lanes 2, 6, 8 and 9. In lanes 3, 5 and 7 no homo/heteroduplexes are 
present. The photograph of the ethidium bromide-stained gels is 
displayed as reversed image. 


Table 2 
Genetic diversity between amplicons subjected to HMA. The length and 
the percentages of genetic similarity between homoduplexes are shown. 


Identification Sequence Similarity (%) FCoV 
13LVI_04a 320 bp 317/320 (99%) Type I 
13LVI_04b 348 bp Extra 5’ 28nt 

19LVI_04a 320 bp 320/320 (100%) 

19LVI_04b 348 bp Extra 5’ 28nt 

22LVI_04a 347 bp 344/347 (99%) 

22LVI_04b 347 bp 

25LVI_04a 320 bp 290/320 (90%) 

25LVI_04b 320 bp 

1OLVI_05a 348 bp 316/320 (98%) 

1OLVI_O5b 320 bp Extra 5’ 28nt 

11LVI_05a 348 bp 303/320 (94%) 

11LVI_O5b 320 bp Extra 5’ 28nt 

33LVI_05a 320 bp 319/320 (99%) 

33LVI_05b 348 bp Extra 5’ 28nt 

99LVI_05a 320 bp 319/320 (99%) 

99LVI_O5b 376 bp Extra 5’ 56nt 

1FLVI_O06a 320 bp 299/320 (93%) 

1FLVI_O6b 320 bp 

32LVI_06a 170 bp 168/170 (98%) Type II 
32LVI_O6b 170 bp 


household were included in branch A indicating the 
presence of different viral strains within each confined 
cat population and within the same individual host. 
Sequences 83 LVI_05, 88LVI_05, 9OLVI_05, 94LVI_05 and 
96LVI_05 came from cats sharing the same household 
than the source animal of sequence 65LVI_05 [74LVI_05, 
81LVI_05, 86LVI_05, 87LVI_O5]. However sequence 
65LVI_0O5 was positioned outside group A, in a separate 
branch with 80% of bootstrap value. The above-men- 
tioned sequences were obtained from a multi-cat 
household where a clinical FIP episode (65LVI_0O5) had 
been reported. Although all the other animals were FCoV 
positive, none developed FIP. Sequence 11LVI_O5b was 
also included in branch A but sequence 11LVI_05a, 
obtained from the same cat was positioned in branch B. 
Sequence 25LVI_04a was grouped in branch A though it 
was obtained from the same animal as sequence 
25LVI_04b placed in a different branch. These two viral 
quasispecies (25LVI_04a/b) came from the same house- 
hold of the source cats of sequences 22LVI_04a/b and 
23LVI_04, which in this analysis were found distributed 
in outlier branches supported by 100% of bootstrap 
value. Sequences 20LVI_O6Z (stool), 2OLVI_O6PBMC 
(blood) and 2OLVI_O6L (ascites) amplified from the same 
animal, were placed in two separate branches, one with 
100% bootstrap value (20LVI_LO6PBMC and 20LVI_O6L). 

Type II FCoV phylogenetic tree included previously 
assigned FCoV Type II sequences and two Portuguese 
sequences obtained during this study, one positive to viral 
quasispecies (32LVI_06; 32LVI_06a/b). All the sequences 
clustered in one branch, although FIPV 79-1146 and 79- 
1683 grouped together with 100% of bootstrap value. Due 
to the low number of nucleotides considered in the align- 
ment (170b), a likelihood mapping analysis was preformed 
(TREE-PUZZLE) to assess its phylogenetic content, revealing 
62.86% of fully resolved, 0% of partially resolved and 37.14% 
of unresolved quartets demonstrating therefore its reliable 
use to estimate phylogenetic relationships. 


4. Discussion 


Until now, no previous study on the presence and 
prevalence of FCoV Types I and II in domestic cats had 
been done in Portugal. In our analysis 47.5% of the 
samples were positive in the RT-PCR assay (57/120) and 
were obtained either from healthy animals or from cats 
with FIP. Moreover, the vast majority of the amplicons 
was Classified as Type I. Type II was found only in two 
samples both from cats with FIP. The values found for 
FCoV Types I and II were in accordance with the 
distribution reported in Europe and USA (Addie et al., 
2003; Benetka et al., 2004) although in Japan and Taiwan 
a higher prevalence of FCoV Type II was observed in 
animals with FIP (Hohdatsu et al., 1992; Lin et al., 2009). 
In our study even though the two FCoV Type II sequences 
were obtained from sick cats, the low incidence of this 
viral biotype does not allow any relationship between 
FCoV II and symptomatic disease. 

The evaluation of FCoV genetic diversity revealed 21% 
of the amplicons positive for viral quasispecies. These 
viral variants were found within the same animal but also 
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Fig. 2. Phylogenetic analysis of FCoV Types I and II sequences based on the 3’-end of the S gene. The trees were constructed by Neighbor Joining based on the 
genetic distances between the nucleotide sequences, calculated by Kimura two-parameter using a transition/transversion rate of 3.23 for Type I sequences 
and 9.99 for Type II sequences. A boostrap analysis of 1000 replicates was performed and values above 80% are indicated on the branches. Branches A-C 
include FCoV sequences with the lowest value of nucleotide substitutions within each branch. Sequences marked with an asterisk (*) were amplified from 


animals with FIP. 


in samples collected from animals sharing a multi-cat 
environment. This finding supports the viral internal 
mutation theory (Dye and Siddell, 2007) and confirms the 
animal exposure to divergent viral strains. The signifi- 
cance of FCoV quasispecies had already been studied 
(Battilani et al., 2003) suggesting a correlation between 
high viral genetic diversity within the host and develop- 
ment of disease. However from our data it was not 
possible to correlate viral diversity with virus subtype or 
disease therefore a more extensive study should be 
performed to assess this connection. 

The phylogenetic analysis of FCoV Types I and II 
sequences revealed an assorted distribution of sequences 
in Type I tree, showing high genetic diversity of the FCoV 
Portuguese sequences, its divergence toward European 
and Japanese sequences (Fq1_1099UK, H9_1100UK, and 
strains KU2, Black and UCD1) confirming the previous 
findings. 

Interestingly, FCoV Type II phylogenetic tree although 
with a limited number of Portuguese and previously 
characterized Type II sequences, showed a homogeneous 


distribution of sequences as also reported by Lin et al. 
(2009) in a similar study, indicative of high similarity 
between distantly related sequences (UK; Netherlands; 
Portugal). The genetic diversity of the Portuguese FCoV 
sequences confirmed and strengthens the evidence that 
this virus persists within the host as mingled viral 
populations whose interaction with the animal may or 
not result in severe disease. 

In conclusion this study provided information on the 
FCoV molecular epidemiology in Southern Europe which 
was unavailable until now and revealed high genetic 
diversity between the Portuguese and previously char- 
acterized FCoV sequences. 
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